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ABSTRACT 

We present a causal Hamiltonian quantum mechanics in 2n- dimensional phase space 
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1. Quantum Contextuality . The most important difference between quantum mechanics and 
a classical stochastic theory is that quantum probabilities are inherently and irreducibly con- 
text (i.e., experimental arrangement) dependent. For any complete commuting set (CCS) of 
observables A, the quantum state \ip) specifies the probability of observing the eigenvalues a 
as | (ct\ip) | 2 , if A were to be observed . If B is another CCS with eigenvalues (3, but [A, B] ^ 0, 
the analogous probabilities \{{3\^)\ 2 , if B were to be measured refer to a different context or 
experimental situation. Each context corresponds to the experimental arrangement to mea- 
sure one CCS of observables. Due to this inherent context dependence quantum mechanics 
does not specify joint probabilities of noncommuting observables. Moreover, the context 
dependence is irreducible, i.e., quantum mechanics cannot be embedded in a classical con- 
text independent stochastic theory. This is the lesson from decades of work, e.g. Gleason's 
theorem 1 , Kochen-Specker theorem 2 , Bell's theorem 3 (where the contextuality corresponds 
to violation of "local realism") and Martin- Roy theorem 4 (which is a direct phase space 
proof relevant to the present work). The contextuality theorems circumscribe the extent 
to which dynamical variables in quantum mechanics can be ascribed simultaneous 'Reality' 
independent of observations. 

2. De Broglie-Bohm . The De Broglie-Bohm (dBB) causal quantum mechanics 5 has shaped a 
paradigm of causal quantum mechanics in which the Position Observable occupies a favoured 
status of a "beable" 6 with values independent of context or observation, whereas other ob- 
servables may have context dependent values. The state of the individual system is char- 
acterized by {\ip(t)),x(t)} where \ip(t)) and x(t) are the state vector and the configuration 
space coordinates, whereas an ensemble of these states corresponds to \ip(t)). For a many 
particle system with the quantum Hamiltonian 

the individual Xi{t) move according to 

dx ' — * 

(Pi)dBB = m i—jjT = ^iS(x(t),t), 

where = Rexp(iS/h), with R and S being real and x denoting (x\, X2, ■ ■ •)• This means 
that the phase space dynamics is determined by the causal Hamiltonian 

and corresponds to the phase space density 

PdBB(x,p,t) = \^(x,t)\ 2 5(p - PdBB (#,*))• 

(We shall set h — 1 henceforth). Integration over momentum shows that the ensemble posi- 
tion density agrees with \if)(x, t)\ 2 . However, as pointed out by Takabayasi 7 , integration over 
position does not yield the quantum momentum density \ip(p, t)\ 2 , where ip is the Fourier 
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transform of ip. This disagreement exhibits the context dependence of momentum in dBB 
thoery: the preexisting momentum probability density given by the dBB theory is assumed 
to be converted into the correct quantum density by means of a measurement interaction 
appropriate to the context of a momentum measurement. On the other hand, position mea- 
surements simply reveal the existing position. Thus the position measurement interaction 
does not play the same role of altering the existing probability distribution. 

3. Motivations For A Causal Quantum Mechanics More Realistic than De Broglie-Bohm 
Theory . The asymmetrical treatment of position and momentum constitutes the breaking 
of a fundamental symmetry of quantum theory which has sometimes been considered as a 
defect of the dBB theory (Holland, Ref. 5, p. 21). We recently constructed a causal quantum 
mechanics 8 ' 9 in one dimension, in which Takabayasi's objection as well as the asymmetric 
treatment of position and momentum are removed. Without invoking the measurement 
interaction, the new causal theory reproduces both position and momentum probability 
distributions of usual quantum theory, and is therefore more realistic than the dBB theory. 

Can we formulate a notion of a maximally realistic causal mechanics (for spinless particles 
with a configuration space of n dimensions) which respects quantum contextuality theorems? 
A mechanics which yields Hamiltonian evolution of phase space variables with a positive 
definite phase space density will be called a 'Causal Hamiltonian Mechanics' or a 'Causal 
Mechanics' in brief. A 'Causal Mechanics' which simultaneously reproduces the quantum 
probability densities of the maximum number of different (mutually noncommuting) CCS of 
obserables as marginals of the same phase space density will be called a 'Maximally Realistic 
Causal Quantum Mechanics'. The definition is nontrivial because the contextuality theorems 
do not allow probability distributions of all possible CCS to be simultaneously reproduced. 

What constraints can we impose selfconsistently on the phase space probability density 
p(x,p,t) of an ensemble of phase space points of a causal theory? Motivated by the success 
in one dimension 8,9 we may require that the quantum position and momentum probability 
densities are reproduced as 'marginals', i.e., 



where ip denotes the Fourier transform of the wave function if). The probability interpretation 
necessitates the positivity condition, 



which rules out many phase space distribution functions such as the Wigner function. 10 
Moreover, positive distribution functions obtained by smoothing the Wigner function 11 do 
not in general reproduce the correct marginals. Neverthless, as emphasized by Cohen and 
Zaparovanny 12 , the uncertainty principle does not preclude the existence of a phase space 
density obeying conditions (1) - (3). A simple example is 




(1) 




(2) 



p{x,p,t) > 0, 



(3) 



Po (x,p,t) = \iJ(x,t)\ 2 \mt)\ 2 - 
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For a causal theory a further condition is necessary if the phase space density is to arise from 
an underlying Hamiltonian dynamics, viz. the "Liouville condition" , 

i«M + g(*£ + Ag)-* (4) 

where a dot denotes time- derivative. As is well known, this condition is an immediate 
consequence of the phase space continuity eqn. 

^ wM -» <*> 

and the existence of a causal Hamiltonian H c (x,p,t) such that 

dH c /dpi = ±i, dH c /dxi = -pi. (46) 

The achievement of de Broglie and Bohm 5 was to construct a causal mechanics obeying 
(1), (3) and (4). The mechanics we constructed 8 ' 9 for n — 1 is more realistic because it obeys 
Eq. (2) in addition. The new mechanics has the phase space density 

p(x,p,t) = \i;(x,t)\ 2 mp,t)\ 2 6 ( [ P dpW,i)| 2 - r dx'\^{ex\t)f) (5) 

\J — oo J — oo / 

where e = ±1, and we have omitted the tilda denoting Fourier transform (i.e., ip(p, t) actually 
stands for ip(p,t) = (p\tp(t))). Eqs. (1), (2) and (3) are obviously satisfied; further, it also 
has a c-number causal Hamiltonian of the form 8 ' 9 

H c =^-(p-A(x,t)) 2 + V(x,t) 

with two quantum potentials A and V (instead of just one in the de Broglie-Bohm theory) 
which depend on the wave function ip. Hence the Liouville condition (4) is also obeyed. 
We shall see that in higher dimensions, we can and should be even more ambitious. 

4. Maximally Realistic Causal Quantum Mechanics . The purpose of the present work is to 
construct a new mechanics which we tentatively call "maximally realistic causal quantum 
mechanics" in 2n-dimensional phase space. In this theory a single phase space probability 
density reproduces the quantum probability densities of n + 1 different CCS of observables 
in spite of the fact that no two sets are mutually commuting. The pleasant surprise is that 
not only the quantum probability densities of position and momentum, but also those of 
n — 1 other CCS of observables can be simultaneously realized. The choice of the n + 1 
CCS whose probabilities are simultaneously realized is not unique; the appropriate choice 
can depend on the context. Different contexts have the same wave function but different 
phase space probability densities; thus the wave function is not a complete description of 
these probabilities. 

The phase space quantum contextuality theorem of Martin and Roy 4 proves that it is 
impossible to realize quantum probability densities for all possible choices of the CCS of 
observables as marginals of one positive definite phase space density. We conjecture that the 
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simultaneous realization of quantum probability densities of more than n + 1 different CCS 
is impossible and hence that the causal theory presented here is maximally realistic. 

Let the state of the individual system be characterized by {\ip(t) >,x(t)} or {\ip(t) > 
,p(t)} since any n independent phase space coordinates are now on equal footing. Due to 
the freedom of canonical transformations, we may assume without loss of generality that 
the CCS (Xi, ■ ■ ■ X n ) is among the n + 1 CCS whose quantum probability densities are 
reproduced in the new causal theory. We assume (without any fundamental justification) a 
one-to-one relation between coordinates and momenta, as this played a crucial role in our 
construction of a causal hamiltonian in one dimension 8 ' 9 . The phase space density must then 
be of the general form, 

n 

P (x,p,t) = w&ttfnsfa-pjfat)) (6) 

which returns the correct marginal \ip(x, t)\ 2 on integration over the momenta. We shall now 
show that the functions p(x, t) can be chosen so as to reproduce the quantum probability 
densities of a 'chain' of n + 1 different CCS, e.g. 

(Xi, X 2 , ■ ■ ■ , X n ), (Pi, X 2 , • • • , X n ), (Pi, P 2 , X 3 , ■ ■ ■ , X n ), ■ ■ ■ , (Pi, P 2 , ■ ■ ■ , P n ), (7) 

where each CCS in the chain is obtained from the preceding one by replacing one phase space 
variable by its canonical conjugate. In the one dimensional case the rquirement of one-to- 
one relation between x and p(x, t) yields two discrete solutions p(x, t) (corresponding to 
e = ±1 in Eq. (5)) which are non-decreasing and non- increasing functions of x respectively. 
Analogously, in the n-dim case, there is a 2-fold ambiguity in determining the phase space 
variables of each CCS in the chain (7) in terms of the preceding CCS, and hence 2 n discrete 
solutions p(x,t). These solutions can be read off from the (^-functions in the 2 n phase 
space densities, each of which reproduces the desired n + 1 quantum probability densities as 
marginals: 

n n 

p(x,p,t) = nm^t)\ 2 u^Aj). (8) 

i=0 j=l 

Here, each fli denotes phase space variables corresponding to one CCS : 

= (pi,P2, ■ ■ ■ ,Pi,x i+ i, ■ ■ ■ ,x n ), for \<i<n-l, 

fin = (pi,P2,---,Pn), 

and the ip(Qi,t) denote appropriate Fourier transforms of ip(Q ,t). Each S(Aj) serves to 
determine flj in terms of flj-i , 

/Pi reixi 
\i>(p'i,x 2 ,- ■ ■ ,x n ,t)\ 2 dp'i- / \ip(e 1 x[,x 2 , - ■ ■ ,x n ,t)\ dx'i, 
-co J —oo 

/Pi 
\fp(Pl, ■ ■ ■ ,Pj-l,P'j,Xj+l, ■ ■ ■ ,£„,t)| 2 G^ 
-oo 

- 1 3 3 \${p\,---,Pj-i,eirfj,Xj+ 1 ,---,x n ,t) ] ?drf j , for 1 < j < n, 
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/Pn 
\^(PU---,Pn-l,P' n ,t)\ dp' n 
-oo 

\tp(Pl,---,Pn-l,e n x' n ,t)\ dx' n , 

-oo 

with 

€t = ±1, for 1 < i < n. 

Since there are 2 n possible choices of the ei, • • • , e n we have here 2 n phase space densities. 
Direct integration over n variables, (using the n 5- functions), yields 



J p(x,p,t)dtii = Ma,f)| : 



which are the correct marginals. Here f2j denotes the n-tuple of phase space variables 
complementary to i.e., = (xi, x 2 , ■ • ■ , Xi,p i+ i, • ■ • ,p n ), with = (j>i, • • • ,p n ) an d 
&n = (xi, %2, ■ ■ ■ , x n ). The condition A ± — determines pi in terms of xi, x 2 , • ■ • , x n and 
t, i.e., • • • , x n , t); A 2 = determines p 2 in terms of pi and x 2 , ■ ■ ■ , x n , t and hence 

p 2 (xi, ■ ■ ■ x n ,t) after substituting pi = pi, and so on. Hence all the momenta are determined 
via the S(Aj) in terms of x±, ■ ■ ■ ,x n ,t, and the phase space density (8) can be rewritten in 
the form (6). (The (^-functions S(Aj) can of course be used equally well to determine the 
coordinates in terms of momenta. E.g. A n = yields x n in terms of Pi, • • • ,p n ,t; Ai-i — 
yields x n -i in terms of p±, ■ • • , p„_i, x n , t and hence in terms of pi, ■ • • , p„, t after substituting 
for x n , and so on). 

The phase space density (8) corresponds to the choice (7) of the n + 1 CCS. The form (6) 
is however more general since it results for any choice of the chain of n+1 CCS which includes 
(X u X 2 ,---,X n ). E.g. for n = 2, all the three chains {{X U X 2 ), (P u X 2 ), (P 1: P 2 )},{(X 1 ,P 2 ), 
(X 1 ,X 2 ),(P 1 ,X 2 )}, and {(X 1 ,X 2 ),(X 1 ,P 2 ),(P 1 ,P 2 )} will lead to Eq. (6), of course with 
different functions p(x,t). 

Consistency Condition on Velocities due to Schrodinger Eqn. The density p(x, p, t) of the 
ensemble of system points depends on the n+1 marginals \ip(0, i: t)\ 2 . Hence the velocities 
of the individual system points will also be constrained by the time dependent Schrodinger 
Eqn. We work out these constraints explicitly when the Schrodinger Eqn. is of the form : 

ihd^/dt = (^2 7^ + U(xfj i) 

in terms of the chosen coordinates and momenta. Starting from the general form (6) of the 
phase space density and taking a partial derivative w.r.t. t with x and p fixed we obtain, 

—p(x,p,t) = 

\ / J — L 

(9) 




dt 

The time dependent Schrodinger Eqn. yields the probability current conservation eqn. 

<9 n 8 

-|^(f,t)| 2 + E^(^t) = 0, (10) 
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where 



ji(x,t) = Re 



in) 



Further, the conservation of the number of system points in phase space yields the phase 
space continuity eqn. (4a) for dp/dt, with Xi = ii(x,p,t), pi = pi(x,p,t). Substituting Eqs. 
(4a) and (10) into Eq. (9) we obtain, 



E 

i=i 



d 

— (vi\^(x, t)\ 2 -ji(x, £)) 
d J [dpi dpi(x,t) 



H 6 (pj -Pj& t)) 
j j=i 



dpi \ \dt dt J 

where v deotes the system point velocities. Thus 

v(x, t) = x(x,p, t) 



x,t)\ 2 l[5(p 3 -p J (x,t))\ = 0, 

3=1 



(12) 



_d_ 

dx; 



(vi\ij(x,t)\- 



d 



+ E 



p=p(x,t) 

dpk d \ 



(13) 



dxi ^ dxi dp k/ 



(ii\ij>(x, t)\' 



p=p 



and 



dt 



x,t) dpi(x,t) ™dpi(x,t). _ 
+ E^I x k (x,p,t) 



dt 



k=l 



dxi 



p=p 



Integrating Eq. (12) after multiplying by dp 1 ■ ■ ■ dp n we obtain the constraint on velocities, 



Umx,t)\ 2 - 3l (x,t)) = o. 

i=l OX i 

Similarly, multiplying Eq. (12) by Pkdpi ■ ■ ■ dp n and integrating, we obtain 

( dp k dp k (x,ty 



(14) 



dt 



dt 



p=p(x,t) 



0, 



which is identically satisfied. 

For n — 1, Eq. (14) implies that the dBB velocity is the unique solution if we wish to 
avoid singularities of the velocity at nodes of the wave function. For n > 1, Eq. (14) can be 
solved for the velocities to yield, 



I m//- + m2 ^dW u (x,t) 
{Vi-v itB )\4>(x,t)\ — . 



where 



-w a , 



and ViB denotes the dBB velocity 



Vi, B = ji(x,t)/\i/j(x,t)\ : 



(15) 
(16a) 
(16b) 
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The velocities given by (15) differ from the dBB velocities due to the term involving the 
antisymmetric tensor W. Eq. (15) was derived directly from Eq. (10) by Deotto and Ghirardi 
and by Holland in their search for atternatives to dBB trajectories 13 . Our derivation shows 
that the Schrodinger Eqn. places no other constraints, for example on pj(x, t). We shall now 
see that in general W has to be non-zero in order that a causal Hamiltonian exist. 

Partial Differential Equations For Velocities From Existence of the Causal Hamiltonian . In 
addition to the constraint (14) due to the Schrodinger Eqn., the velocities must also obey 
partial differential eqns. which follow from the requirement of existence of a c-number causal 
Hamiltonian H c (x,p,t). If X{ and p\ are derived via Hamilton's equations, 

/ -» ,\ • I ( dH c {x,p,t) \ 
Vk(x,t) = x k \^ = f ^ j _ (17) 



dpi(x,t) dpi dH c (x,p,t) 



v=v dxj 



p=p 



(18) 



dt dt 
Defining 

H c (x,t) = H c (x,p,t)\ p ^ m , (19) 
and substituting Eqs. (17) and (18), we have, 



dH c (x, t) _ ( dH c dp k dH c 



dxi \ dxi dxi dpk 



p=p 



dp t (x,t) + dPu Vk ^ ty (20) 



dt dxi 

In order that a function H c (x,t) obeying the partial differential eqns. (20) exist, the inte- 
grability conditons 

d 2 H c d 2 H c 

(21) 



dxidxj dxjdxi 

must hold. Substituting (20) into (21) we obtain the n(n— l)/2 conditions on the velocities, 

(1 < i < j < n), 

d d d 

where 

d d 

fij & *) = -fapi *) ' -fa.Pi *) • (23) 

Note that these partial differential eqns. for the velocities for existence of a causal Hamil- 
tonian are derived without any assumption about the functional form of the Hamiltonian. 
When we substitute Eqs. (15) for the velocities (given by Schrodinger Eqn.) into Eqs. (22), 
we obtain n(n — l)/2 partial differential eqns. for the n(n — l)/2 functions Wa(x,t): 

d_ ( hi dw M \ _ d_ ( f kl dWjA 

dx t \mx,t)\ 2 dx e ) dx 3 \mx,t)\ 2 dx e j + ij u ' 1 ] 
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where 

dfi d d 

F ij = -fif + Q^(fkjV k ,B) ~ g^(fkiVk, B ). (25) 

Except for very special wave functions (e.g. factorizable wave functions), the pi determined 
to fit quantum probability distributions of n + 1 CCS yield ^ and ^ 0, and Eqs. (24) 
do not have the trivial solution W k e = 0. In contrast, in the dBB theory which reproduces 
only the position probability density (pi)dBB = diS which yields fij = = and hence 
Eq. (24) is obeyed with Wke = 0. Thus, for general wave functions, departure from dBB 
velocities (i.e., Wke ^ 0) is needed for existence of a causal Hamiltonian when we insist on 
reproducing quantum probability distributions of n + 1 CCS (n > 1) of observables. 

Determination of Causal Hamiltonian . With velocities so determined, H c is found by inte- 
grating Eqs. (20) without any assumption about the form of the Hamiltonian H c . We shall 
show that a causal Hamiltonian exists without making any claim of its uniqueness. We now 
make the ansatz, 

H c (x, p, t) = ± ^~^ (f?t))2 + V(x, t). (26) 

i=l % 

We find from Hamilton's eqns. (17), 

Ai(x, t) = pi(x, t) - rriiVi (26a) 
which yield the Ai\ we then calculate V(x, t) from, 

V(x, t) = H c (x, f) - £ (Pi ~ Ai)2 . (266) 

i 1 

This completes the determination of a causal Hamiltonian (26) which contains n+1 quantum 
potentials. The explicit forms of the potentials for n — 1 have been given in Ref. 9. For 
n = 2 they follow from the velocity formula (15), and Eqs. (25) and (28). 

Explicit Formulae For Velocities For n = 2 . In 2 dimensions Eq. (24) simplifies to a first 
order partial differential eqn. 

d dW^ d dW^o 

£-(212)^-^12 = 0, (27) 



OX 2 OX i OX i ox 

where 

9n(x,t) = f 12 (x,t)/\^(x,t)\ 2 . (27a) 

Since no time derivatives occur, t may be considered as a fixed parameter and Eq. (27) 
solved by Lagrange's method by considering a curve x\ = xi(s), x 2 = x 2 (s) with 

dx\ dx 2 dW\ 2 

dg 12 /dx 2 dg X2 /dxi F 12 

On this curve 



^12 n dW_ 



0, 



12 



ds ds 



'12- 
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Inserting back the fixed parameter t, we obtain the most general solution, 



/ 

Jo 



( 



dg 12 /dx 2 



12 



) 



(x 1: x 2} t)d,Xi 



(28) 



where (i) the argument x' 2 along the path of integration equals x 2 at x[ = X\ and is determined 
for other values of x[ by the condition dg\ 2 = along the curve, and (ii) h is an arbitrary 
function of g 12 and t. The velocities given by (15) will contain a degree of arbitrariness 
corresponding to the choice of the function h(gi 2 ,t). 

5. Conclusions . The new causal quantum mechanics in n dim. configuration space has the 
following important properties, (i) It reproduces quantum probability distributions of n + 1 
CCS of observables with one positive definite phase space density, (ii) It has a c-no. causal 
Hamiltonian which contains n + 1 quantum potentials, (iii) It has velocities (and hence 
trajectories) which (for general wave functions) differ from dBB velocities for n > 1 due to 
the insistence on reproducing quantum probabilities of n + 1 CCS. The velocities contain 
some arbitrariness (e.g. the function h in 2 dim.) in spite of these constraints, (iv) It has 
position- momentum correlations in individual events (given by p(x,t)) different from dBB 
theory. Applications to quantum chaos and possible experimental tests need further work. 
We also hope to compare numerically the trajectories implied by the present work with those 
given by the de Broglie-Bohm theory in a future communication. 
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